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ABSTRACT 

We investigate optical and near-IR color gradients in a sample of 172 low-inclination galaxies spanning 
Hubble types SO-Irr. The colors are compared to stellar population synthesis models from which 
luminosity-weighted average ages and metallicities are determined. We explore the effects of different 
underlying star formation histories and additional bursts of star formation. Our results are robust 
in a relative sense under the assumption that our galaxies shared a similar underlying star formation 
history and that no bursts involving more than ~10% of the galaxy mass have occurred in the past 
1-2 Gyr. Because the observed gradients show radial structure, we measure "inner" and "outer" 
disk age and metallicity gradients. Trends in age and metallicity and their gradients are explored 
as a function of Hubble type, rotational velocity, total near-IR galaxy magnitude, central surface 
brightness, and scale length. We find strong correlations in age and metallicity with Hubble-type, 
rotational velocity, total magnitude, and central surface brightness in the sense that earlier-type, 
faster rotating, more luminous, and higher surface brightness galaxies are older and more metal-rich, 
suggesting an early and more rapid star formation history for these galaxies. The increasing trends 
with rotational velocity and total magnitude level off for V Iot > 120 kms -1 and Mjf < —23 mag 
respectively. This effect is stronger for metallicity (than age) which could reflect a threshold potential 
above which all metals are retained and thus metallicity saturates at the yield. Outer disk gradients 
are found to be weaker than the inner gradients as expected for a slower variation of the potential 
and surface brightness in the outer parts. We find that stronger age gradients are associated with 
weaker metallicity gradients. Trends in gradients with galaxy parameters are compared with model 
predictions: these trends do not agree with predictions of semi-analytic models of hierarchical galaxy 
formation, possibly due to the effect of bar-induced radial flows. The observed trends are in agreement 
with chemo-spectrophotometric models of spiral galaxy evolution based on CDM-motivated scaling 
laws but none of the hierarchical merging characteristics, implying a strong dependence of the star 
formation history of spiral galaxies on the galaxy potential and halo spin parameter. 
Subject headings: galaxies: spiral — galaxies: color gradients — galaxies: photometry — galaxies: struc- 
ture 



1. INTRODUCTION 

Color gradients in galaxies reveal information about 
the nature of their stellar populations via age and metal- 
licity trends, and the amount and distribution of dust. 
The technique of using broad-band colors as a probe 
of the stellar populations and star formation histories 
(SFH) of galaxies, pioneered by, e.g. Searle, Sargent, 
& Bagnuolo (1973) and Tinsley & Gunn (1976), is still 
reflected in the modern studies of Peletier & Balcells 
(1994), de Jong (1996), and Bell & de Jong (2000; here- 
after BdJOO). Early analyses, however, suffered sig- 
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nificantly from degeneracies between age and metallic- 
ity. Worthey (1994) quantified the age-metallicity de- 
generacy that exists in optical broad-band colors as 
Aage/AZ ~ 3/2. This implies that the composite spec- 
trum of an old stellar population is nearly indistinguish- 
able from that of a younger but more metal-rich popula- 
tion (and vice versa). This degeneracy can be partially 
broken with infrared photometry (e.g. H or K band) 
in addition to optical colors (de Jong 1996). Cardiel 
et al. (2003) have quantified the relative ability of dif- 
ferent color and absorption-line index combinations to 
constrain physical parameters of composite stellar pop- 
ulations. Their results find that inclusion of an infrared 
band improves the predictive power of the stellar popula- 
tion diagnostics by ~ 30x over using optical colors alone. 
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The interpretation of broad-band color gradients also re- 
lies on a careful mapping of the dust extinction within 
a galaxy (Witt, Thronson, & Capuano 1992; de Jong 
1996). While dust opacity is much reduced at redder 
wavelengths, its effects may still be non-negligible in the 
determination of the stellar content of late-type spiral 
galaxies and must be considered in the final interpreta- 
tion of color gradients. 

A successful galaxy formation theory must determine 
the formation timescales for bulges and disks and repro- 
duce the observed stellar ages and chemical distributions 
in galaxies as a function of galaxy properties. Competing 
bulge/disk formation scenarios entail distinct structural 
and kinematic signatures, star formation rates (SFRs), 
and mixing of the bulge and disk material, that a study 
of age and metallicity indicators could help disentangle. 
Current theories for the formation of disk bulges gen- 
erally fall under two basic pictures: monolithic collapse 
(e.g. Larson 1974; Calberg 1984) and hierarchical merg- 
ing (e.g. Kauffmann 1996; Cole et al. 2000). The mono- 
lithic dissipative collapse of a gas cloud at high redshift is 
consistent with strong bulge metallicity gradients (on the 
order of 0.5 dex/decade within the central ^2 kpc, but 
flattening in the outer parts (Larson 1974)), with metal- 
licity decreasing with increasing radius, shallow positive 
age gradients (the centers are slightly younger than the 
outer parts), and old stellar populations. In this "mono- 
lithic" scenario the disk would accrete onto the already 
formed bulge, and thus be younger than the bulge. The 
dominant signatures in the hierarchical scenario depend 
on the time evolution of the merger rate: if mergers were 
most frequent at high redshift, as predicted by present- 
epoch cosmological models, signatures should be similar 
to those of dissipative collapse. More recent mergers and 
tidal disruptions would result in flatter metallicity gradi- 
ents, and a significant fraction of young or intermediate- 
age stars. In either of the above two scenarios, secular 
evolution can also drive disk gas and stars to the cen- 
ter via bar instabilities. In general, this scenario pre- 
dicts that bulges will be similar in metallicity and age 
and structurally linked to disks (Courteau, dc Jong, & 
Broeils 1996; Courteau 1996). Note, however, that the 
number of bar formation - gas inflow - star formation 
- bar dissolution cycles that can occur in a galaxy may 
be limited, as the central mass formed in the first cycle 
can prevent subsequent formation of a bar Wyse (1999). 
1 Radial migration of disk material (stars and gas) can 
also be induced by non-axisymmetric spiral waves (Sell- 
wood & Binney 2002), which would necessarilly flatten 
metallicity gradients in disks. 

Current observational results (O/H gradients, SFRs, 
constant IMF, gas distribution) favor an inside-out for- 
mation for the Milky Way, but the details about age 
and metallicity gradients in our own Galaxy are con- 
flicting (Freeman & Bland-Hawthorn 2002). Integrated 
studies of stellar populations in external galaxies can re- 
strict the range of possible interpretations, but forma- 
tion timescales and chemical evolution of external spiral 
galaxies have barely been explored (e.g. Matteucci 2002). 

1 Only progressively larger bars in the centers of exponential 
bulges would be allowed to form in a recurring scenario as a result 
of the disrupting dynamical effect of a growing nucleus (H.-W. Rix 
1999, private communication, as reported in Carollo 1999). 



These issues have become even more important now that 
cosmologically-motivated numerical simulations are ca- 
pable of resolving the formation of disks and bulges. 
Furthermore, semi-analytic galaxy evolution models use 
recipes to describe bulge and disk evolution, but if dif- 
ferent types of galaxies have different bulge formation 
mechanisms, these models may give misleading results 
(e.g. Somcrvillc & Primack 1999; van den Bosch 2000), 
unless a variety of bulge formation mechanisms are mod- 
eled (e.g. Cole et al. 2000). 

Existing studies of spiral bulges and disks have placed 
tentative constraints on the source of their color gradi- 
ents and which galaxy parameters most strongly corre- 
late with SFHs and, hence, on specific formation scenar- 
ios, but some of the results are conflicting. In an study 
of color gradients in early-type bulges (<Sb) using HST 
and ground-based data, Peletier et al. (1999) infered that 
early-type bulges are old (absolute ages not well deter- 
mined) with a small age spread (< 2 Gyr) for all early- 
type bulges measured at the bulge effective radius. The 
old and narrow age range makes it unlikely that these 
bulges would have formed via secular evolution. The 
same conclusion does not hold for their (few) later-type 
galaxies. Peletier et al. (1999) conclude that the intrin- 
sic color gradients of early-type bulges are caused mainly 
by metallicity gradients (in agreement with other simi- 
lar studies; e.g. Mchlcrt et al. 2003), consistent with 
the monolithic collapse scenario. However, based on the 
lack of r 1 ' 4 shaped bulges found in a later study of bulge 
SB profiles of SO-Sbc galaxies using HST near-IR high 
resolution imaging, Balcells et al. (2003) have revised 
their interpretation and conclude that these bulges could 
not have formed from violent relaxation in mergers, but 
rather are more consistent with secular formation pro- 
cesses. On the other hand, in an HST study of spatially 
resolved colors of high redshift (z ~ 0.5) galaxies in the 
Hubble Deep Field, Abraham et al. (1999) rule out metal- 
licity gradients as major contributors to galaxy color dis- 
tributions. They do find, in agreement with Peletier et 
al. (1999), that large bulges are significantly older than 
their disks and therefore rule out secular evolution for- 
mation processes in favor of a gradual disk formation by 
accretion of gas. However, the flat metallicity distribu- 
tion of the bulge is not easily explained in this scenario 
(particularly since the lack of "r 1 / 4 " bulges rules out ma- 
jor merging as a cause of gradient flattening). 

The most comprehensive study to date of color gradi- 
ents in disk galaxies is that of BdJOO, who compute age 
and metallicity gradients for a sample of 121 nearby low- 
inclination SO-Sd galaxies. They conclude that galaxy 
color gradients are due in most part to gradients in age 
and metallicity in their stellar populations (in the sense 
that inner regions are older and more metal rich than 
outer regions) and contend that dust reddening mainly 
affects the metallicity gradients, but is likely too small 
to affect their conclusions. In comparing the SFHs with 
galaxy parameters, they conclude that the SFH of a 
galaxy is primarily driven by surface density and that 
the total stellar mass of a galaxy is a less important pa- 
rameter that correlates significantly with metallicity, but 
not age. Kauffmann et al. (2003a) also conclude that the 
recent SFHs, as probed by the Usa absorption line in- 
dex and the 4000A break, of over 100,000 galaxies are 
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more strongly correlated with surface mass density than 
stellar mass. BdJOO argued that these correlations could 
be the result of a surface density-dependent star forma- 
tion law, coupled with galaxy mass-dependent chemically 
enriched gas outflows. Bell & Bower (2000) further ex- 
plored this idea by constructing a family of simple models 
for spiral galaxy evolution for comparison with the obser- 
vational trends in SFH with galaxy parameters. Indeed 
they found that the data are consistent with the propo- 
sition that the SFH of a region within a galaxy depends 
primarily on the local surface density of the gas but that 
additional ingredients, such as galaxy mass dependent 
mfall(outfiow) of primordial(metal-enriched) gas and/or 
formation epoch, are required to fully explain the obser- 
vational results. 

In spite of genuine progress in recent years in the study 
of spectro-photometric gradients in spiral galaxies, our 
ability to model and interpret them in terms of forma- 
tion models is still limited due to the lack of data and 
the intricacies involved in converting observed quantities 
to reliable ages and metallicities. In this study, we use 
the largest catalog of deep optical and NIR galaxy col- 
ors to date to revisit the comparison of broad-band color 
gradients with stellar population models using a range 
of SFHs and basic assumptions about the dust distri- 
bution. We follow the approach developed by BdJOO, 
exploring additional parameter ranges and using an ex- 
tended data base by combining the BdJOO data with our 
own collection of deep, optical and NIR surface bright- 
ness profiles (Courteau, Holtzman, & MacArthur 2004, 
in prep.; hereafter Paper III), for a total sample of 172 
galaxies. For the combined data base we determine lo- 
cal average ages and metallicities in radial bins for each 
galaxy and compute gradients in age and metallicity as 
a function of radius. We pay special attention to the 
effects of fitting out to different physical extents for in- 
dividual galaxies and the distinction between inner and 
outer galaxy gradients. In particular, we find that false 
trends can be inferred if the radial extent of the gradi- 
ent fits is not taken into account, which in turn leads to 
erroneous conclusions about the galaxy parameters that 
drive their SFHs. 

The outline of the paper is as follows. In |J21 we de- 
scribe the data base from which color gradients are com- 
puted. Our color gradients are presented in Paper III, 
and further details of the BdJOO sample can be found in 
their §2. In Jllwe explore the range of galaxy parame- 
ters in our sample and their intrinsic correlations which 
must be considered when comparing trends in age and 
metallicity gradients with galaxy parameters. In jj^l we 
discuss the stellar population models to be compared to 
the data and the different star formation prescriptions 
adopted. Dust models and its potential effects on our re- 
sults are discussed in Optical-NIR color-color profiles 
and their matching population models are presented in 
f|51and trends with galaxy parameters are explored. The 
technique by which the data are fitted to the SSP models 
is presented in jjZI Local and global age and metallicities 
are p rese nted in H8.1l and their gradients are discussed in 
M8.3I and H8.4I and contrasted with previous results from 
BdJOO. We conclude with a discussion of the mechanisms 
that control stellar evolution in spiral galaxies and com- 
pare our results with existing models of galaxy evolution 
in P 



2. THE DATA 

Our data base from which color gradients are mea- 
sured is a combination of the data in BdJOO and the 
compilation by Courteau et al.(2004, in prep., Paper III) 
of 1063 digital images in the B, V, R, and H pass- 
bands of 324 nearby late-type spiral galaxies collected 
at the Lowell Observatory and Kitt Peak National Ob- 
servatory. For the current analysis we consider only the 
face-on and moderately tilted galaxies (i < 60°) in the 
sample. All galaxies were selected from the Uppsala Gen- 
eral Catalogue of Galaxies (UGC, Nilson 1973) to have: 
Zwicky magnitude vtlb < 15.5, blue Galactic extinction 
A B = 4 x E(B-V) < 0f5 03), and blue major axis 
< 2'.2. For the computation of homogeneous structural 
parameters and color gradients we use the isophotal map 
determined at i?-band and applied onto all other images 
(BVH) of a galaxy. The SB profiles are reliable (with SB 
errors < 0.12 mag arcsec -2 ) down to ~ 26 mag arcsec -2 
for optical bands and ~ 22 mag arcsec -2 at if-band. 
For the purpose of this color-based analysis, we further 
require that the galaxies have measured surface bright- 
nesses in at least two optical and one near-IR (H) band 
observations out to at least 1.5 7J-band disk scale lengths 
and with SB errors of less than 0.12 mag arcsec -2 . This 
leaves us with 51 galaxies (25 type I, 17 type II, 9 tran- 
sition) 2 with extended reliable color gradients. 

Radial profiles were corrected for Galactic extinction 
using the reddening values of Schlegel et al. (1998). We 
do not attempt to correct our SB profiles for internal 
extinction, but discuss the possible effects of dust on our 
results in [JS] The SB profiles were degraded to the worst 
seeing FWHM before computing color gradients. For a 
full description of the data see Paper III. 

In order to increase the signal-to- noise (S/N) of our 
color profiles, the SB profiles were averaged into radial 
bins scaled by the NIR-band disk scale length, /inir, and 
we required at least 3 data points per bin. The bin sizes 
were defined as follows: 
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The second, finer, grid (grid2) was used to see if the 
coarser binning of gridi hides any important features in 
the color profiles. The difference in using the two differ- 
ent binnings is small. In particular, this test confirmed 
that the measured gradients do not depend on the central 
pixels, which are likely to be more affected by a central 
concentration of dust in spiral galaxies. Hence we sub- 
sequently use gridi only (as this is the binning of the 
BdJOO sample). 

2 For type I SB profiles (Freeman 1970), the inner profile always 
lies above the SB of the inward extrapolation of the outer disk, 
whereas type II systems have a portion of their brightness profiles 
lying below the inward disk extrapolation. We define a transition 
case for luminosity profiles that change from type II at optical 
wavelengths to type I in the infrared. See Paper I for a definition 
of the different profile types. 
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Table 1. 



Hubble Type 


Sa 


Sab 


Sb 


Sbc 


Sc 


Scd 


Sd 


mean H — K 


0.30 


0.26 


0.30 


0.28 


0.28 


0.24 


0.21 



Note. — Mean 2MASS H - K isophotal colors for 
non-barred galaxies. Adapted from Figure 16 of Jarrett 
et al. (2003). 



In order to look for trends in color gradients with 
galaxy parameters, we collected, for as many galaxies as 
possible, their morphological type, absolute magnitude, 
ii-band disk central surface brightness (CSB), ii-band 
disk scale length, inclination, total ii-band magnitude, 
and rotational velocity I^ot (see below). The structural 
parameters (disk scale length, h, and CSB, (Iq) were com- 
puted using our bulge-to-disk decomposition scheme (see 
Paper I). Total galaxy magnitudes encompass the light 
out to SB levels given above, and include an extrapo- 
lation of the SB profile to infinity with an exponential 
fit to the last quarter of the SB profile 3 . The extrapo- 
lated magnitude increment (0.07 iJ-mag on average) is 
added to the isophotal magnitude which is corrected for 
Galactic extinction (Schlegel et al. 1998). Stellar masses 
are computed using the prescription of Bell & de Jong 
(2001). 

We merged our sample with the 121 galaxies from 
BdJOO which spans a wider range in Hubble type (SO-Irr) 
and includes a few low surface brightness (LSB) galaxies. 
For the computation of color gradients, BdJOO required 
at least two optical (BVRI) and one if-band observa- 
tions per galaxy. The data reductions and radial binning 
(using only grid\) are very similar to our own. For 3 
galaxies in common with both data sets, there is excellent 
agreement in the overall calibration and radial profiles at 
least for the BVRH bands (Paper III). 

In order to compare the two data sets directly when 
looking for trends in age and metallicity with galaxy pa- 
rameters, we converted our ii-band magnitudes and sur- 
face brightnesses into iC-band using the 2MASS H — K 
color transformations derived from Jarrett et al. (2003). 
Since we are interested mainly in qualitative trends in age 
and metallicity with galaxy parameters, coarse transfor- 
mations are adequate and we took the mean H—K values 
for each Hubble type, derived from their Figure 16 and 
listed here in Tabled The modified H magnitude will 
be denoted as H m , to represent the iT-band equivalent of 
the ii-band measurement. We do, however, consider the 
disk scale lengths in the H and K to be directly compara- 
ble. Mollenhoff & Heidt (2001) find hj/h K = 1.01±0.19 
for a sample of 40 spiral galaxies and the relationship will 
be even tighter for fin/hx, hence a direct comparison of 
the H and K band disk scale lengths is justified. 

It has been suggested that age and metallicity gradi- 
ents may depend upon the potential of the galaxy and 
thus vary as a function of V^ot (Dalcanton & Bernstein 
2002; Garnett 2002a). In order to search for such a trend 
we augmented out catalog with H I or Ha line widths 

3 This accounts for outer disk truncations that would not be fit 
properly using the entire profile from our B/D decompositions in 
Paper I. 



(which we take to be W50/2, half the line width at 50% 
peak flux) collected from the literature. These data will 
enable us to test for a distinct signature in the color pro- 
files from the more prominent dust lanes seen in galaxies 
with Vrot > 120 kms" 1 (Dalcanton & Bernstein 2002) 
and/or if there is a flattening of the metallicity of galax- 
ies above V TQ t ~ 150 kms -1 (Garnett 2002a). The line 
widths were collected from various sources and transfor- 
mations were derived to ensure uniformity between them 
(see Appendix ^ . Differences between measurements 
for a given galaxy for different sources were typically less 
than 10%, which is sufficient for our purposes. 

In order to ensure sufficiently accurate V TO t measure- 
ments, we restricted our sample to galaxies with incli- 
nation greater than 35°. We are left with 98 galaxies 
in our sample that have reliable rotational velocity mea- 
surements. 

3. GALAXY PARAMETERS 

Trends in age and metallicity (local values and gradi- 
ents) with galaxy parameters, may bear an imprint from 
int rinsi c correlations within the sample itself. In Fig- 
urc lAll we plot the galaxy parameters used in this analy- 
sis (h, no, MH m ,K, and V ro t) against each other. The 
most notable correlations are the luminosity-rotation 
speed relation, Mh^^k versus log(Kot) (the so-called 
"Tully-Fisher" relation), and the size-luminosity rela- 
tion, Mh^^k versus \og(h). All correlati ons with are 
weaker and have large scatter. In Figure lA~2l we plot the 
4 basic galaxy parameters as a function of morpholog- 
ical type. Trends are seen with galaxy parameters for 
T > 4 such that later-type galaxies have smaller rota- 
tional velocities and total magnitudes, and lower central 
surface brightness. On the other hand, there is no clear 
trend of disk scale length with Hubble type (and hence 
the trend with total magnitude is largely driven by the 
trend with CSB, Courteau et al. 2003). Data points are 
few for T < 4, but deviations from the T > 4 trends 
seem to occur. A larger spread is (barely) noticeable 
in V I0 t, but is further supported by the larger spread in 
Mh^^k for which we have more galaxies (unfortunately, 
reliable Vr 0t s could not be found for the earliest-type 
galaxies.) The increasing trend in /io with decreasing 
T seems to level off at T < 4, but the spread is still sig- 
nificant. Finally, earliest- type galaxies have the shortest 
scale lengths, as expected. We shall return to these cor- 
relations in the analysis of trends in age and metallicity 
gradients (^Hl- 

4. STELLAR POPULATION MODELS 

In order to translate color information into constraints 
on the underlying ages and abundances of the stellar 
populations, the color gradients must be compared with 
stellar population synthesis models. In their most basic 
form, commonly referred to as simple stellar populations 
(SSPs), these models provide evolutionary information 
for a coeval population of stars born with a given compo- 
sition and initial mass function (IMF). The closest physi- 
cal analog to such SSPs are globular cluster systems from 
which the SSP models are calibrated (e.g. Schiavon et al. 
2002). Several such SSP models have been produced by 
a number of independent groups and are in a constant 
state of modification as improvements to many of the in- 
put parameters (e.g. stellar libraries, model atmospheres, 
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convection, mass loss, mixing) come to light. There are 
discrepancies among the different models that, depend- 
ing on the application, may result in significantly differ- 
ent interpretations of observations of unresolved stellar 
populations. In order to determine if these model dis- 
crepancies could affect our analysis, we have compared 
two independent sets of SSPs; the 2003 implementation 
of the Bruzual & Chariot (2003) models (hereafter re- 
ferred to as GALAXEV), and the Project d'Etude des 
Galaxies par Synthese Evolutive (PEGASE 2.0) models 
of Fioc & Rocca-Volmerange (1997). 

The stellar populations of interest here are those of spi- 
ral galaxies spanning the full Hubble sequence (SO-Irr). 
The approximation of a single evolved stellar population 
clearly does not apply to evolved, complex stellar popu- 
lations of spiral galaxies. However, if one assumes that 
galaxies are composed of a superposition of SSPs, born 
at different epochs, rates, and metallicities, one can use 
the SSPs to develop model grids that mimic the range 
of plausible galactic stellar populations and these can be 
compared with galaxy color profiles. The stellar popula- 
tion model grids are created by taking the single burst 
SSPs, with constant stellar IMF and fixed metallicity, 
and convolving them with a given SFH. A few simpli- 
fying assumptions are inherent in this formulation. We 
are assuming that the IMF does not change as a func- 
tion of time or galactic environment. The validity of this 
assumption may be questioned from a theoretical per- 
spective as the IMF is expected to vary systematically 
with star formation environments (Chabricr 2003, and 
references therein). However, significant observational 
evidence of star formation in small molecular clouds, rich 
and dense massive star-clusters, as well as ancient metal- 
poor stellar populations, reveals a remarkable uniformity 
of the IMF (Kroupa 2002). The use of fixed- metallicity 
SSP implies no chemical evolution. This is clearly an 
unrealistic assumption with regards to galaxy evolution 
(as it also ignores feedback effects from stellar winds and 
supcrnovae), but it still allows us a reasonable compar- 
ison of relative metallicities and ages within and among 
galaxies (Abraham et al. 1999; BdJOO; Gavazzi et al. 
2002). 

The specific choice of IMF, parameterized as 
£(logm) oc mT x , can also potentially affect our results. 
To gauge the importance of this effect on our color-based 
analysis we compared model grids obtained with some 
of the most widely used IMF characterizations; (a) the 
single-slope Salpeter (1955) with 

a; =1.35, 0.1 < m/M Q < 100, (2) 
(b) the Kroupa (2001) IMF (PEGASE models only) 
with: 

£0 = 0.3, 0.1 < m/M Q < 0.5 
£1 = 1.3, 0.5 < m/MQ < 120, (3) 
and (c) the Chabrier (2003) single-star Galactic disk 

IMF (GALAXEV models only) which is parameterized 

as: 



£(logm) oc 



exp 



(log m— log m c ) 
2^ 



-1.3 



0.1 < m/M G < 1; 

1 < to/M q < 100 
(4) 



The Chabrier (2003) IMF is preferred over the Kroupa 
(2001) and Salpeter (1955) IMFs because of its better 
agreement with number counts of brown dwarfs in the 
galactic disk and its theoretical motivation. 

The integrated spectrum, F\(t), of a stellar population 
with SFR *f?(t) is computed as the convolution of a single 
burst stellar population of given metallicity, f\(t), with 
the given SFR: 



F x (t)= / *(t-t')f x (t')dt'. 
Jo 



(5) 



The left panel in Figure lA~3l compares the model tracks 
from the GALAXEV models with the Salpeter (1955) 
(eq. 0) and Chabrier (2003) (eq. 0) IMFs, and the right 
panel compares those from the PEGASE models for 
Salpeter (1955) and Kroupa (2001) (eq. EH) IMFs. The 
difference between the various IMFs on the model grids 
are small, thus the specific choice of IMF (within current 
observational constraints) will not affect our results. We 
opt to use the Salpeter (1955) IMF for the remainder 
of this analysis (primarily to facilitate comparison with 
other studies). 

Also over-plotted in the right panel of Figure IA3I are 
the GALAXEV model with Salpeter IMF tracks for di- 
rect comparison with the PEGASE model tracks for the 
same IMF (although note that the upper mass cut-off is 
100 M in the GALAXEV models and 125 M in the PE- 
GASE models). While differences in the grids between 
the different models are noticeable (most likely due to 
different treatments of the thermally pulsing asymptotic 
giant branch phase, see e.g. Maraston 2003), they are 
only significant in the bluest regions where the SB errors 
of our observations are also large (i.e. the outer regions of 
our spiral galaxies where sky subtraction errors become 
large, are, for the most part, bluer). Slight differences in 
the age and metallicity gradients would be inferred us- 
ing the different models, but we are mainly interested in 
relative quantities, and thus are not concerned with the 
small absolute differences between models. For the rest 
of this analysis we refer only to the GALAXEV models. 

We have explored two different SFR regimes, both pa- 
rameterized by a star formation timescale r. The first is 
the simple exponential SFR, ^ cxp (t), 



xp(0 — ^ 
T 



-t/l 



where 



l- e -A/r- 



(0 



(7) 



The second is the so-called "Sandage" SFR, ^'san(i) 
first parameterized by Sandage (1985), 



*Sa„(<) = c s 



and 
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where c oxp and cg an are normalization constants (to 1 



M Q ). Figure IA4I shows the time evolution of the expo- 
nential (top panel) and Sandage SFRs (bottom panel). 
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The exponential SFR falls off at a given rate for all pos- 
itive values of r. However, in order to cover the full 
color-color space spanned by the observations, negative 
values of r were also included which allow for increasing 
SFRs. The Sandage SFR is characterized by a delayed 
rise in the SFR, followed by an exponential decline, the 
rate of which is determined by the value of r. For values 
of r > A, where A is the age of the galaxy (fixed to 13 
Gyr in this work), the SFR is still rising at the present 
time. The average age of a stellar population of a given 
t is computed as: 



(A)=A 



(10) 



For the exponential SFR, this equates to, 



1 - e" A / T fl + A/t) 
(A) exp = A-T—, [ _ ± A/ T) (11) 
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Figurc lASl shows the resulting model grids for the expo- 
nential (blue) , and Sandage (red) SFRs overlaid on each 
other. Lines of constant average age (A) are the dashed, 
roughly vertical lines (labeled in Gyr), while lines of con- 
stant metallicity are solid and run closer to the horizontal 
(labeled in Z, the mass fraction in elements heavier than 
helium, where Z = 0.02). The overall shape of the two 
model grids are quite similar, the main difference being 
that the Sandage SFH has a younger average age for a 
given B — R and the metallicity is slightly higher for a 
given R — H. The Sandage grids also cover a wider range 
in ages (which would increase the magnitude of our mea- 
sured age gradients). 

The reported radial age gradients inferred from color 
gradients (e.g. BdJOO) seem to be at odds with the lack 
of age gradient observed in open clusters of the Galaxy 
(Freeman & Bland-Hawthorn 2002). Would it be pos- 
sible that the apparent radial age gradients observed in 
spiral galaxies are simply due to small amounts of young 
frostings of star formation on top of an underlying old 
population? In order to test for such effects in our age 
determinations, we produced new model grids for the 
same underlying SFH but with an additional burst of 
star formation at different times and fractional masses. 
The burst was added as a Sandage (1985) type profile 
(eq. |EJ with r = 0.1 4 . The average age for the models 
with a burst superimposed on the underlying SFR is now 
computed as, 

(A) = A - a-^/o^^W^ + ^/^^d* (13) 

(i-b f )j A ^(t)dt + b f j t U b (t')dt 

where bf is the mass fraction of the burst, tb is the 
time since the burst, i.e. the burst of star formation oc- 
curs tb Gyr before the current age of the galaxy, A, and 

4 This timescale is justified on the basis that star bursts are con- 
sistent with a constant SF lasting 10-100 Myr (Meurer 2000). 



t' = t — Figures lA"6l and VK1\ demonstrate the effect of 
such a burst of star formation superimposed on an un- 
derlying exponential SFR with ages (tb) of 1 and 4 Gyr 
respectively. The left panel in both figures shows the ef- 
fect if 10% of the galaxy's total mass was involved in the 
burst, and the right panel shows a 50% (by mass) burst. 
The models certainly do not cover the same extent in the 
color-color space for recent bursts. With only 10% of the 
total galaxy mass in the burst (Fig. IA6I left panel) the 
grid does not extend to the red as far as those without 
bursts. The effect is strongest (~ 0.3 mag) in the optical 
colors. Also, the burst grid does not extend as far into the 
blue optical colors (< 0.1 mag) for the most metal-rich 
galaxies. This is because we are essentially adding more 
older stars to the SFHs that are rising at 13 Gyr (r < 
for V&ezp, or r > 13 Gyr for ^>san', This effect is stronger 
at higher metallicity because the evolution of B — R color 
with time steepens for higher metallicity stellar popula- 
tions at ages ~ 1-3 Gyr in the GALAXEV models.) The 
lines of constant age also steepen. "Frosting" also shifts 
the R — H colors to redder values (~0.1 mag) at a given 
metallicity. This, in fact, prevents the burst grid from 
reaching red-ward enough as required by the observa- 
tional datafor the inner-most parts of the galaxies (see 
Figs. IA9HA13|) . Either the bursts are accompanied by 
significant amounts of dust, or 1 Gyr old bursts can be 
ruled out for the central par ts of late- type galaxies. The 
50% by mass burst in Figure lA"6l right panel, shows that 
the grid covers only a narrow range in color-color space 
that does not extend far enough into the red or the blue 
to agree with the observations. For older bursts of star 
formation, the model grids have colors fully consis tent 
with the no-burst models. This is shown in Figure E3 
for a pure exponential SFH (blue) overlaid with an expo- 
nential plus tburst = 4 Gyr Sandage-like burst (red) . The 
grids are almost identical with respect to their shape and 
location (with slight narrowing of the 50% mass burst 
grid). The only difference is the average ages (eqs. EH 
& I13|) , being younger for a given color of the burst grid. 
Thus, adding a burst of star formation that is older than 
~ 1 Gyr has a similar effect as changing the overall form 
of the SFH; the grid shape and color coverage are the 
same, but the average age can change significantly from 
one scenario to another. Not only do the age determina- 
tions change, but the age gradients can also be different 
for a given SFH (co mpa re Sandage and exponential SFH 
model grids in Fig. IA5|I . Since the actual SFHs of real 
galaxies are unknown, the absolute values of the gradi- 
ents measured with this color-based technique cannot be 
trusted. We can, however, trust that a gradient exists in 
a galaxy, but not its magnitude. We can also compare 
relative gradients among different galaxies if we assume 
self-similar SFHs and no significant bursts of star forma- 
tion within the last 1-2 Gyr. 

A reasonable estimate of the validity of this assump- 
tion can be obtained from the study of Kauffmann et 
al. (2003a) who used the HSa absorption line index and 
the 4000A break age indicators along with the Bruzual 
& Chariot SSPs to constrain the SFHs, dust attenua- 
tion, and stellar masses of over 100,000 galaxies from the 
Sloan Digital Sky Survey (SDSS). For each galaxy they 
estimated Fburst, the fraction of the total stellar mass 
formed in bursts in the past 1-2 Gyr. Their Figure 5 
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shows the fraction of galaxies with F(, urs i > at the 50% 
and 97.5% confidence levels in stellar mass bins ranging 
from 8.0-12.0 logio(M ). From this and estimates of the 
galaxy stellar masses (Bell & de Jong 2001), we infer the 
fraction of galaxies in our sample that have undergone 
bursts of SF in the past 1-2 Gyr is - 1% (~ 12%) at 
the 97.5% (50%) confidence level. The fraction of galax- 
ies that have Fb urst 3> 10% will be even smaller and we 
conclude that our assumption of no recent burst is valid. 

5. DUST EFFECTS 

While the presence of dust in late-type spiral galaxies 
is well-established, its distribution and effective optical 
depth remain poorly constrained. Extinction by dust 
could potentially mimic a color gradient and its effects 
must be considered in determining the stellar content of 
galaxies from colors. We attempt to quantify the effects 
of dust extinction on colors with simple dust geometry 
extinction models, adopting the MW extinction curve 
and albedo values from Gordon, Calzetti, & Witt (1997). 
The extinction vector for the simplistic foreground screen 
model with Ay = 0.3 is show n (as the arrow) in the up- 
per left corner in Figures IA9HA131 The vector is in the 
same general direction as the observed gradients of the 
galaxies, but given the unrealistic nature of this model, 
we do not consider it further (but include it for visual 
com parison) . Also shown in the top left corner in Fig- 
ures IA9IIA131 is the more realistic "face-on triplex" dust 
model of Disney, Davies, & Phillipps (1989) and Evans 
(1994). This model assumes that the stars and dust have 
exponential distributions in both the radial and vertical 
directions with radial scale lengths ft* and hd, and ver- 
tical scale lengths z* and z&- Although scattering is not 
taken into account, Byun, Freeman, & Kylafis (1994) and 
de Jong (1996) have shown that its effects are negligible 
in face-on galaxies (as photons are just as likely to be 
scattered into the line of sight as they are out of it.) We 
therefore use the dust absorption curve to compute the 
dust optical depth; 

r A (0) = ^ry(0)(l-a A ) (14) 

Ty 

where t\ and a\ are the dust optical depth and albedo 
respectively. 

Studies of multi-wavelength data that solve simulta- 
neously for the intrinsic color of stellar populations and 
reddening and extinction effects from dust find central 
optical depths, ry (0), in the range 0.3-2.5 for types Sab- 
Sc (Peletier et al. 1995; Kuchinski et al. 1998; Xilouris 
et al. 1999). The HST and ground-based color gradient 
study of early-type bulges (<Sb) by Peletier et al. (1999) 
reveals dust extinctions of Ay = 0.6-1.0 mag (or Ah = 
0.1-0.2 mag) in galaxy centers but negligible extinction 
beyond one effective radius. Multi-band SB profile mod- 
eling of massive edge-on disks suggests that the dust is 
confined to a thin extended plane, with Zd/z* < 0.7 and 
h d /K ~ 1.4 (Xilouris et al. 1999; Matthews & Wood 
2001). Masters, Giovanelli, & Haynes (2003) studied in- 
ternal extinction in the near-IR for a sample of 15,244 
2MASS galaxies by looking at the inclination dependence 
of various photometric parameters, and concluded that 
galaxies with M A - 51og(h) > -20, -20.7, and -20.9 in 
the J, H, and if s -bands respectively are spared any ex- 
tinction, but that disk opacity increases monotonically 



with disk luminosity. Application of the triplex model 
for these galaxies favors small dust-to-star scale height 
ratios (zd/z* ~ 0.5, in good agreement with Xilouris et 
al. 1999) and face-on central opacities of t°(0) = 0.7,0.3 
at H and K s respectively. 

We adopt the same vertical and radial scale length 
ratios for the dust and stars; using Zd/z* — 0.7 and 
hd/h* = 1.4 instead stretches the gradients by only 0.018 
in B - R, 0.019 in R - H, and 0.021 in R - K, and the 
central F-band optical depth to Ty(0) = 1 (or pole-to- 
pole V-band optical depth of 2). The higher end of the 
measured range of Ty( 0) — 2 is shown in the upper left 
corner of Figures lAltil & IaTTI (right panels) . 

The triplex models lie parallel to, and thus could 
contribute significantly to, the observed gradients (note 
there is no overall calibration in the dust models so they 
can slide to any region on the plot). For the gradients to 
be entirely due to dust requires extremely high central 
optical depths (tv(0) ~ 5) to reproduce the data. Also, 
the triplex models alone cannot reproduce the significant 
color gradients observed from the half-light radius (de- 
noted by open circles) outward. Thus, while dust is likely 
a contributor to color gradients, age and metallicity ef- 
fects must also be invoked. 

Dalcanton & Bernstein (2002) have analysed optical 
and infrared color maps of 47 extremely late-type (bulge- 
less) edge-on spirals spanning a wide mass range (40 < 
V c < 250 kms -1 ). They find that higher mass (rota- 
tion velocity) galaxies have more prominent dust lanes 
and have redder colors than the stellar population grids 
(whereas the less massive galaxies have colors consistent 
with the population grids). Analysis of their R—K s color 
maps, and a comparison of colors between their edge-on 
sample and the face-on sample of de Jong (1996) suggests 
that dust plays little role in all but the most massive 
galaxies in their sample (V c > 120 kms -1 ). 

Previous analyses have ruled out dust effects on the 
color gradients in galaxies based on the fact that they 
see no correlation between the gradients and inclination. 
Dust effects on gradient profiles as a function of bulge-to- 
total ratio (B/T) and inclination were modeled by Byun, 
Freeman, & Kylafis (1994). Their Figure 6 shows B - 
I versus radius for model galaxies with different B/T 
ratios and Ty (0) ranging from 0.0-10.0 for inclinations 0- 
85°. In the 0-50° range the differences in the profiles are 
very small and would not be detected in the observations. 
Our measured radial gradients also show no inclination 
dependence, but we do not consider this grounds to rule 
out significant dust extinction. 

Based on optical-IR imaging of SO-Sbc galaxies, 
Peletier & Balcells (1996) found that "dust-free" col- 
ors of galaxy disks are not significantly different from 
their bulges. They derived bulge colors from minor- 
axis wedges in images of early- type edge-on spirals. The 
wedges are presumed dust-free above the disk plane. 
They conclude that (dust-free) bulge and disk colors 
are very similar with A(B — R) = 0.045 ± 0.097 and 
A(R - K) = 0.078 ± 0.165. We find a significantly dif- 
ferent result; as can be seen in Figure ^8] (compare with 
Fig. 2 in Peletier & Balcells 1996), our bulges are much 
redder than their disks with A(B - R) = 0.29 ± 0.17 
as A(R - K) = 0.30 ± 0.17. This could be due to gen- 
uine dust extinction in our galaxies or there is a funda- 
mental difference between our respective samples and/or 
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analysis methods. Their disk SB profiles were measured 
along 10°-wide wedge apertures centered 15° away from 
the disk major axis, to avoid the prominent dust lanes 
near the major axis of their inclined galaxies. Naturally, 
this technique is sensitive to a vertical disk color gradi- 
ent. Dalcanton & Bernstein (2002) have suggested that 
all thin disks are embedded in a red stellar envelope. 
Whether this envelope is redder or bluer than the thin 
disk depends on the rotational velocity of the galaxy: 
redder envelopes for V ro t < 120 kms -1 and bluer oth- 
erwise, but the color of the envelope is similar for all 
disk galaxies. The redder thin disks of the galaxies with 
V rot > 120 kms -1 are attributed to strong dust lanes 
observed in their B — R color maps, which disappear in 
the V ro t < 120 kms -1 galaxies. Clearly, these red stellar 
envelopes, and hence the presence of vertical color gradi- 
ents in disk galaxies renders the interpretation of Peletier 
& Balcells (1996) "disk" colors difficult. The 10° offset 
from the major axis may result in a measurement of the 
"red envelope" stars which are intrinsically older than 
the thin disk stars (though not necessarily redder since 
thin disk stars could be reddened by a central dust lane) . 

While optical-IR color imaging alone is not sufficient 
to break the degeneracy between dust and stellar pop- 
ulation effects on the color gradients, the tentative con- 
sensus to date is that dust is generally not a significant 
contributor to galaxy colors in low-mass/low-luminosity 
spiral galaxies, but is likely important in more mas- 
sive/brighter galaxies. It must be borne iin mind that 
future studies of high-resolution IR and FIR imaging 
and absorption-line spectroscopy may radically alter this 
view. 

6. COLOR-COLOR PROFILES 

Optical-NIR color-color profiles are show n fo r the 
Cour teau et al. sample of in Figures fA9HA10l and IA12T 
IA13I (see Figs. 1-5 in BdJOO for similar plots with their 
sample and our Fig. IA11I which is only a subset of the 
BdJOO sample with available V TO t valuess). Typical ob- 
servation errors a re shown as crosses in the lower right 
corners of Figures IX"9TlA13l From left to right the crosses 
represent calibration errors, and average sky subtraction 
errors for the innermost and outermost bins. The galaxy 
centers are indicated by solid symbols: circles for type-I 
galaxies, triangles for type-II, and asterisks for transi- 
tion galaxies (as defined in Paper I). Open circles denote 
the half-light radius of the disk. The line types for the 
galaxy color profiles represent bar strength: solid for bar- 
less (A), dashed for mild bars (AB), and dot-dashed for 
strong bars (B). It is conceivable that the type-II "dip" 
could be due to dust extinction and/or stellar popula- 
tion effects, possibly linked to the presence of a bar, and 
or inner-disk truncation (see Paper I). However, we do 
not see any distinction between the different SB profile 
types in the color gradients. This argues against dust 
as a major factor for the type-II phenomenon as dustier 
galaxies would be redder and have larger gradients. No 
conclusion can be drawn about the effects of a bar due to 
the small size of our sample and the fact that non-barred 
galaxies may have once harbored a bar that dissolved 
after mixing the stellar pop ulation. 

Over-plotted on Figures IA9HA131 are our convolved 
Bruzual & Chariot (GALAXEV) stellar population mod- 
els with an exponential SFH at different metallicities. 



Lines of constant average age (A) (see eq. I10|) are the 
dashed, roughly vertical lines, while lines of constant 
metallicity are solid and nearly horizontal. In the upper 
left corner, triplex and foreground screen dust models 
are plotted with the galaxy center and disk half-light ra- 
dius denoted by solid and open circles respectively. All of 
the galaxies show significant color gradients that are con- 
sistent with gradients in both age and metallicity with 
the central parts ha ving o ld er me an ages and being more 
metal rich. Figures IA10I & IA11I for the Courteau et al. 
and BdJOO samples reveal a largely consistent picture of 
the radially-resolved colors of spiral galaxies, though the 
Courteau et al. galaxies extend to higher metallicity (0.1- 
0.2 mag in R— H). The origin of this discrepancy is not 
clear, but we note that the optical versus J — K colors of 
galaxies are not well fit by single metallicity models (e.g., 
Fig. 15 of Bell et al. 2003). It may well be that, due to 
TP-AGB prescriptions, systematic errors of 0.2 mag be- 
tween H and K predicted colors are inevitable. Thermal 
emission from the telescope at AT-band may also affect 
sky flats of the type used by BdJOO. Bearing in mind 
these potential sources of systematic uncertainty, we feel 
that the degree of similarity between the Courteau et al. 
and BdJOO galaxy samples is satisfying. 

The small number of earlier-type galaxies (Sab-Sbc) 
in our sample makes it difficult to infer any tr ends in age 
and metallicity with galaxy type, but Figure 1X281 shows 
hints that later-types are slightly less metal-rich than the 
earlier types. 

In Figiircs lXlOl fc lATTl we divide the Courteau et al. and 
BdJOO samples, respectively, at V ro t 5* 120 kms -1 , the 
threshold above which all edge-on galaxies have promi- 
nent central dust lanes (Dalcanton & Bernstein 2002). 
Due to inclination restrictions, we could only retrieve re- 
liable H I line- widths for 28 of our 51 galaxies and 70 out 
of the 121 BdJO O samp le. A trend with V ro t is detected 
in Figures IA10I & IA11I with faster rotators being more 
metal rich and having older mean ages; also consistent 
with them having a higher dust content. However, ac- 
cording to the triplex models, the true signature of dust 
is not simply a redder color, but also a str etche d color 
profile (compare triplex profiles in Figures IA10I & IA11I 
for TZ-band optical depths of ry(0) = 1, left p anels , and 
7y(0) = 2, right panels). For our sample (Fig. lAlOjl . the 
length of the gradients appears to be about the same in 
both velocity bins. For the BdJOO sample fFig, IAlT|) the 
gradients could be more extended in the larger V TO t bin, 
but the stretching is not predominantly along the R — K 
axis, as in the triplex dust models. Clearly we cannot 
attribute with absolute certainty the redder colors of the 
faster rotators to dust effects. 



7. MODEL FITTING 

Ages and metallicities are determined by fitting the 
SPS models to the radial galaxy colors using a maximum- 
likelihood approach similar to that of BdJOO (see their §3 
for further details). We compute a finely spaced grid in r 
using equationsEl&El and interpolate (linearly) between 
the SPS metallicities. Treating the colors and errors of 
each annulus separately, we compute an average age and 
metallicity per annulus by minimizing the x 2 statistic: 
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where N is the number of passbands (at least 3; 2 
optical plus 1 IR), (5//tot,i is the total error in passband i, 
and fi c is the best normalization between the model and 
data computed as a weighted average: 
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Errors for the individual age and mctallicity measure- 
ments as well as their gradients are estimated using a 
Monte Carlo approach. One hundred realizations of the 
model fits are obtained for each radial bin using errors 
drawn from a normal distribution of the observational 
errors (which include calibration, sky subtraction, and 
flat fielding errors added in quadrature). For each real- 
ization, gradients are computed as weighted linear fits to 
the parameter (age and metallicity) determinations as a 
function of the radial scale length. The weights in the fits 
are taken as the A\ 2 — 1 interval for each annulus- model 
fit. The error for the measured gradients and individual 
ages and metallicities are taken as half the interval con- 
taining 68% of the 100 Monte Carlo realizations (i.e. the 
la confidence interval). 

A few galaxies in both samples fall outside the model 
grids and thus cannot be fit reliably. These were re- 
moved from the sample for the model fits, but are wor- 
thy of some discussion. As can be seen in Figure lAlH five 
of our galaxies lie significantly red-ward of the model 
grids in their central R — H color. The dust vector in- 
dicates that such red colors could result from a central 
dust concentration. Another explanation could be an ex- 
tremely (unrealistically) metal-rich central stellar popu- 
lation. Three of our galaxies lie considerably red-ward 
of the model grids in their B — R colors. Given the 
current WMAP measurement of the age of the universe 
(13.7 Gyr; Spergel et al. 2003), it would be unrealistic for 
these galaxies to have extremely old central regions. Dust 
likely contributes to the red colors, but they could also re- 
sult from erroneous seeing measurements. In the BdJOO 
sampl e 3 ver y late-type (Sdm/Irr) galaxies (see left panel 
of Fig. lAllI) and the extended tails of another 3 galaxies 
lie blue-ward of the model grids in their R — K colors. 
These low surface brightness galaxies may suffer from cal- 
ibration andsky uncertainties couple with the weakness 
of the model predictions at such low metallicities. Fi- 
nally, the BdJOO sample has 3 early-type galaxies whose 
entire profiles lie red-ward of the grids in their B — R 
colors. Again, extremely old ages are improbable, but 
significant dust absorption could contribute to the red 
colors of these three galaxies. Thus, unlike some of the 
studies cited in |3 these observations remind us of our 
fragile understanding of dust effects and the challenge we 
face in trying to separate them from stellar population 
effects. 

8. RESULTS 

8.1. Local and Global Trends in Age and Metallicity 



Figure IA14I shows the local average age (left) and 
metallicity (right) as a function of local surface bright- 
ness. There is a clear trend in the local age and metal- 
licity as a function of local SB in the sense that regions 
of higher SB are older and more metal rich, but with 
large scatter, as also found in BdJOO. This correlation 
suggests that the local potential plays an important role 
in the S FH. 

Figure IA15I shows an effective average age, (A) eg, for 
a galaxy, taken as the measured mean age in the 0.5- 
1.5ft. averaged radial bin, as a function of central sur- 
face brightness (left plot) and total galaxy magnitude 
(right plot). There is a clear trend in (A) c g as a func- 
tion of both ^o,ff m ,K an d Mjj^,k with more luminous 
and highe r CSB galaxies having older {A) e g (as expected 
from Fig, IA14|l . However, in both cases, the roughly lin- 
ear increasing trend at ^o t H m ,K 18.5 mag arcsec -2 
and MH m ,K ^ —22.5 mag seem to flatten for brighter 
values. For Mn m ,K —22.5 there is large scatter and a 
much weaker (if any) trend of increasing age with total 
magnitude (but still containing the oldest galaxies). In 
Figure lA"l6l we plot (A} e g as a function of hu,K (left) and 
V rot (right). As expe cted from the correlation of Mh^^k 
with Vtat (Fig. lAlfl . a trend of increasing (A) e g with 
V m t is detected with a similar turnover in the slope of 
the correlation at higher rotational velocities (V TO t ^120 
kms -1 , the loc ation of the vertical dotted line). The left 
panel of Figure lAl6l mav show a weak trend of increasing 
age with ft, however, the scatter is very large. Note also 
that the spread in age is larger for smaller ft galaxies, 
and reduces with increasing ft. 

Figures IA17I & IA18I similar trends for the effective 
metallicity, \og(Z c g/Z & ) versus fJ>o t H m ,Ki Mff m ,ic, and 
Vrot, but no trends are seen with ft. Effective metallic- 
ity increases with total magnitude up to MH m ,K ~ ~22, 
but likely saturates for brighter galaxies at roughly solar 
metallicity (with large scatter). Similarly, in the right 
panel of Figure IAT81 \og(Z e g/Z Q ) increases with V rot up 
to ~ 120 kms -1 and then levels off with smaller scatter. 
There is no dependence of \og(Z c g/ Zq) with ft. 

8.2. Global Trends in Age and Metallicity Gradients 

One of the conclusions from BdJOO is that the ampli- 
tude of the age gradients increases from HSB to LSB 
galaxies. However, this is likely an artifact of their lin- 
ear gradient fitting technique out to different number of 
radial bins (scaled by the galaxy disk scale length), any- 
where from 2-4 bins (i.e. out to 1.5-3.5 scale lengths). 
Assuming that the LSB galaxy photometry generally 
does not extend to as many scale lengths as the HSB 
photometry, non-linear gradients in the galaxies (e.g. 
steeper in the inner regions, flattening at large radii), 
could mimic a trend with SB. To demonstrate this, in 
Figures lA 1 91 fc lA"20l we plot average age, (A), and metal- 
licity, \og(Z/Zo), respectively, as a function of radius. In 
each case, the left panels give radius in terms of disk scale 
lengths, and the right panels have the physical radius in 
kpc. Modulo fairly large uncertainties, we see that gradi- 
ents are not linear over the extent of the galaxy. The in- 
ner and outer slopes are often quite different, sometimes 
even changing sign. Also, the galaxies are not all mea- 
sured out to the same number of scale lengths. When 
plotted against number of scale lengths, most galaxies 
have steeper (negative) inner gradients that flatten off, 
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or even turn around, in the outer parts. With the excep- 
tion of extremely early or late-types, the outer slopes are 
quite similar for most galaxies. However, when plotted 
as a function of physical scale (kpc), the slopes can be 
quite different, with a rough trend for the earliest types 
being ste epest, a nd lev eling off for later-types. 

Figures IA19I fc lA"20l demonstrate that fitting gradients 
out to a different number of scale lengths yields mislead- 
ing results. In general, the larger the fit baseline, the 
shallower the measured gradient. We verify this asser- 
tion in Figure IA21I where we plot age gradients fit out 
to different numbers of radial bins for all the galaxies as 
a function of CSB (0-1.5hff,K (top left), 0-2.51ih,k (top 
right), 0-3.51ih,k (bottom left), and 0-4. 5Hh,k (bottom 
right)). Indeed, fitting fewer bins results in systemati- 
cally larger (usually negative, but occasionally positive) 
age gradients. The more extended the baseline, the flat- 
ter the gradient. LSB galaxy profiles do not extend as 
far as HSBs (in scale lengths, as indicated by the lower 
CSB galaxies dropping out of the longer baseline fits) 
and thus appear to have larger gradients. 

Given the dangers of fitting gradients out to differ- 
ent radial extents, in the remainder of this analysis we 
consider "inner" and "outer" gradients fitted over the 
\ st -2 nd radial bins and 2 nd -4* /l bins, respectively 5 . This 
restriction greatly reduces the available galaxy param- 
eter range, especially for the LSB galaxies. Also note 
that inner gradients will be much more affected by dust 
(if present) than outer gradients. 

8.3. Age gradients 

In Figure fS22l we plot the inner (left panels) and outer 
(right panels) gradients measured in disk scale lengths 
(upper panels) and kpc (lower panels) as a function of 
CSB (left plot) and total magnitude (right plot). Unlike 
BdJOO, we do not see a trend of d(A)/dr (Gyr/h) (in- 
ner or outer) with CSB. Note that the low CSB galax- 
ies are missing in the outer gradient plots so one can- 
not assess any trend for the outer gradients with great 
confidence. The current data suggest that outer gra- 
dients are generally smaller than those within. When 
age gradients are plotted as a function of kpc, the in- 
ner gradients are slightly steeper for higher CSB galax- 
ies {/io,H nl ,K 17.5). The 3-4 outliers at high positive 
gradients are the SO galaxies whose color gradients are 
inverted (bluer inward). Such anomalous gradients have 
been seen before in SO galaxies (e.g. Emsellem et al. 2002) 
and are often interpreted as the result of a recent gaseous 
accretion followed by (central) star formation. These are 
also the few galaxies that deviate to much younger ages 
from the general trend of decr easing effective age with 
galaxy type (see left plot of Fig. lA27)l which would agree 
with the above interpretation of recent accretion and SF. 

There is a slight trend (with large scatter) for more 
luminous galaxies to have steeper d(A)/dr (Gyi/h) neg- 
ative gradients (inner and outer) , and the fainter galaxies 
extending to positive gradients. However, this trend dis- 
appears when the gradients are plotted as a function of 
kpc. In this scale the outer gradients are greatest (neg- 
ative or positive) for the fainter galaxies, and become 

5 Note that these two definitions cannot be interpreted as bulge 
and disk gradients as the inner fit goes out to 1.5 disk scale lengths, 
and bulge scale lengths are typically ~0.13 (Paper I). 



negligible for the brighter ones. A similar, somewhat 
stronger trend of steepening gradient (in scale len gths) 
with increasing disk scale length is seen in Figure IA23I 
(left). The trend in d(A)/dr (Gyr/kpc) is reversed (ig- 
noring the 3-4 steep positive gradient points) such that 
bigger galaxies have smaller gradients. This may reflect 
the larger galaxies as being less likely to have had any 
recent activity/interactions, as the small gradients com- 
bined with the older ages for the bigger galaxies suggest 
an overall old age for the largest galaxies. Similar trends 
are seen with V m t (as expected from the correlations be- 
tween Vrot with h and M in Fig. lAljl . but they are weaker 
(largely because of the smaller sample). 

8.4. Metallicity gradients 

Figures IA24I fc 1X251 recast the information presented 
in Figures IA22I fc lA"23l discussed above, but now in terms 
of metallicity gradients. 

Metallicity gradients are less well measured and are 
more sensitive to dust effects than age gradients (accord- 
ing to the triplex dust models discussed in fJSJ) . Neverthe- 
less, the trends in the metallicity gradients with galaxy 
parameters are very similar to those with age gradients 
(see S l8.3f) . The trends in the magnitudes of the gradients 
as a function of galaxy parameters are comparable except 
the correlation changes sign. This suggests a direct anti- 
correlation between age and metallicity gradients within 
a galaxy which is seen in Figure lA"26l (although with some 
scatter), such that stronger age gradients are associated 
with weaker metallicity gradients (and vice versa). 

9. DISCUSSION 
Our main findings can be summarized as follows: 

(i) Our relative age determinations are robust under 

the assumption that the underlying SFR of all disk 
galaxies is similar, and that no major star bursts 
have occurred in the past 1-2 Gyr (which is justi- 
fied on the basis of Kauffmann et al. (2003)). 

(ii) Dust cannot be ruled out as a contributor to color 
gradients in spiral galaxies, although it is unlikely 
that the gradients are largely due to dust extinc- 
tion. 

(iii) Age and metallicity correlate strongly with local 
surface brightness: higher SB regions tend to be 
older and more metal-rich (see Fig. IA14]) . This in- 
dicates that the local potential plays a significant 
role in the SFH of spiral galaxies. 

(iv) Age and metallicity, measured at an effective ra- 
dius of 1ft,, increase with earlier Hubble type, 
Mn m k, Koti and /io k but the trends flatten 
for f < 4, M Hm>K '< -22.5 mag, V rot > 120 
kms " 1 , and h q. h.k < 18.5 mag arcsec -2 (see 
Figs. lA15HA18l fc lA2Tj) . Age also correlates weakly 
with scale length (with scatter decr easing with 
galaxy size, see rig ht plot of Fig. EH - 

(v) Age and, to a lesser extent, metallicity gradients 
show radial structure, with generally steeper gra- 
dients in the inner parts of the galaxy. Care must 
thus be taken when defining the gradient fit region 
(see Figs.EIUMB- 



Structure of Disk Dominated Galaxies: II. Color Gradients 



11 



(vi) Age gradients, as measured as a function of scale 
length, show correlations with l uminosity, size, and 
rotational velocity (Fig. IA22|) . especially in the 
inner disk. Trends in metallicity gradients with 
galaxy parameters are weaker (Figs. IA24I fc IA25JI . 

(vii) Age gradients do not correlate with CSB, contrary 
to the findings of BdJOO, which we attribute to 
inco nsiste ncies in their fit radii (see left plot of 
Fig. IA22)l . This is not at odds with statement (hi) 
since and h are not strongly correlated. 

These observations are consistent with a picture where 
(i) higher surface brightness regions of galaxies formed 
their stars earlier than lower surface brightness regions, 
or at a similar epoch but on shorter timescales, and (ii) 
the SFHs at a given SB level, which lead to age gradients, 
are modulated by the overall potential of the galaxy such 
that brighter/higher rotational velocity galaxies formed 
earlier. These trends reach saturation for the brightest 
and highest CSB galaxies. 

An earlier formation time for more massive galaxies 
is clearly in conflict with hierarchical galaxy formation 
which predicts that more massive galaxies form late. 
This discrepancy could be remedied if there is a mecha- 
nism at work that prevents the gas in lower mass galaxies 
from cooling and forming stars at early times. Feedback 
processes are often invoked as a possible solution, but 
no prescription has been found thus far that agrees with 
all observational constraints (e.g. Kauffmann et al. 2003; 
Bell et al. 2003). 

In a related manner, the semi-analytic models of hier- 
archical galaxy formation of Kauffman (1996) predict a 
correlation between bulge-to-disk ratio and luminosity- 
weighted age, such that the bulges of late-type spirals 
should be older (up to ~ 4 Gyr) than those of early- 
type spirals, although this correlation could be erased by 
any significant inflow from the disk after the last major 
merger. Figure IA27I (left) shows a strong correlation of 
effective age with galaxy type, but in the opposite sense 
of Kauffmann's prediction; later-type galaxies have sys- 
tematically younger effective ages. Clearly, mixing by 
a bar could erase the predicted age trend. If bars trig- 
ger radial flows the gradients of strongly barred galax- 
ies would be flat (the flattening would occur on short 
enough timescales). To look for such mixing effects, we 
plot age and metallicity gradients as a function of barred- 
ness in Figure IA29I No trends are seen when the gradi- 
ents are plotted in scale lengths, but against kpc, the 
strongly barred galaxies may have smaller age gradients 
(inner and outer), though this observation is based on 
small statistics. This agrees with Martin & Roy (1994) 
who found similar flattening of O/H metallicity gradi- 
ents with bar strength (their trend also disappears when 
the gradients are plotted in scale lengths). Thus, under 
all assumptions, our observations cannot be reproduced 
by the Kauffmann (1996) models. The same model defi- 
ciency responsible for the backwards age-size correlation, 
as compared against observations, could also be the cause 
of the discrepancy in bulge morphology versus age seen 
in the Kauffmann (1996) models. 

Prantzos & Boissier (2000; hereafter PB00) present 
chemo-spectrophotometric models of spiral galaxy evolu- 
tion. They adopt a Schmidt-type law for the SFR, which 
is proportional to the gas surface density and varies with 



galactocentric radius (dynamical time), assume an in- 
flow rate (of unenriched gas) that decreases exponentially 
with time and increases with surface density and galaxy 
mass, and assume that the gas settles into an exponen- 
tial disk (bulges are not modeled), but do not consider 
radial flows (e.g. by viscosity which could create gradi- 
ents, or due to a bar which could flatten radial gradi- 
ents on small timescales). The star formation efficiency 
and infall rate are free parameters tuned to match the 
Milky Way observational constraints. The models are 
then extended to other disk galaxies by adopting the 
Cold Dark Matter(CDM)-based scaling relations of Mo, 
Mao, & White (1998). As such, their galaxy disk radial 
profiles are fully described by just two parameters: ro- 
tational velocity, V ro t (assuming a constant disk-to-halo 
mass ratio), and the halo spin parameter, A. A third 
parameter describing the formation redshift would be re- 
quired for a description fully consistent with the CDM 
hierarchical models of galaxy formation but, as of yet, 
there is no clear definition for the time of formation for 
an individual galaxy. Hence, PB00 assume that all galax- 
ies started forming at the same epoch (13.5 Gyr ago) but 
evolve at different rates. 

The PB00 models give predictions for O /H abundances 
and gradients that as could be measured in bright H II 
regions in nearby galaxies. Note that O/H determina- 
tions in H II regions probe the present-day ISM abun- 
dances and are insensitive to abundance evolution with 
time. Still, their model predictions and observations of 
O/H gradients can be compared, at least indirectly, to 
our stellar luminosity weighted gradients. PB00 find that 
the absolute central abundance for a given A > 0.03 cor- 
relates with V TO t (and total magnitude) such that faster 
rotators have larger central abundances, but the trend 
saturates for V xo t ^ 220 kms -1 . Furthermore, at a given 
Vrot, the central abundance decreases with increasing A 
but again saturates for galaxies with V TO t > 220 kms -1 
above which central abundances are high, regardless of 
A. By extension, we should see a trend of increasing cen- 
tral abundance with V TO t, with significant spread (larger 
at lower V I0 t) due to different values of A, but this trend 
would flatten and show less dispersion abov e V m t > 220 
kms -1 . This is roughly what Figures TA17I fc IA18I show 
for the effective metallicity at the half-light radius ver- 
sus yiu m ,K and V ro t (central values are too uncertain 
due to likely higher concentrations of dust and seeing 
mismatches, and PBOO's models do not consider bulges). 
Unfortunately, our galaxies do not exceed V TO t ~ 250 
kms , but the general trend above, including the de- 
crease in scatter with V I0 t, is confi rmed. 

The right side of Figure IA24I shows our measured 
metallicity gradients as a function of total galaxy mag- 
nitude; gradients in disk scale lengths are shown in the 
upper panels while gradients in kpc are shown under- 
neath. The left and right panels show the "inner" and 
"outer" gradients respectively. When plotted against 
scale length, the inner or outer gradients show no clear 
trends. However, when plotted in kpc, a trend emerges 
such that brighter galaxies have smaller gradients. Since 
scale length increases with luminosity (see Fig. lAlfl . the 
correlation goes away when the gradients are measured 
per disk scale length. If this effect is real, the disap- 
pearance of a trend of metallicity gradients in units of 
scale length may suggest a self-similar pattern in disk 
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galaxies. As suggested by Combes (1998), a "universal" 
slope per disk scale length might be explained by the 
viscous disk models of Lin & Pringle (1987), although 
model predictions are not conclusive as of yet. The Gar- 
nett et al. (1997) compilation of O/H gradients shows a 
similar signature, where abundance gradients in dex/kpc 
are steeper and exhibit greater scatter for lower luminos- 
ity disks, but this trend goes away when gradients are 
expressed in dex/h. Note, however, that the abundance 
gradients compiled in van Zee et al. (1998) and reported 
in PBOO do show a clear increase with magnitude when 
plotted as dex//i, thus thwarting the interpretation of a 
universal abundance gradient per scale length. 



We have measured compelling trends for the effective 
age and metallicity in spiral galaxies and their inner and 
outer disk gradients as a function of surface brightness, 
luminosity, rotational velocity, and size. These trends are 
not well reproduced by current semi-analytical models of 
galaxy evolution but will undoubtedly serve as effective 
constraints for future models. An important limitation 
inherent in our color-based analysis is its sensitivity to 
dust extinction which could mimic gradients in age and, 
particularly, metallicity. A spectroscopic analysis of ra- 
dially resolved line indices of spiral galaxies is underway 
in order to alleviate this problem. 



APPENDIX 
H I LINE WIDTHS 

H I line widths were collected from various sources in the literature; Courteau (1997) [7gals]; Bottinelli et al. (1990) 
[74 gals]; Theureau et al. (1998) [11 gals]; Haynes et al. (1999) [1 gal]; de Blok & Bosma (2002) [1 gal]; de Blok, 
McGaugh, & Rubin (2001) [4 gals]. For the Bottinelli et al. (1990) and Theureau et al. (1998) samples, the line width 
at 20% and 50% of the peak flux, W50 and W20, were provided, though not for all galaxies. In cases where only 
W20 was available, a conversion to W50 was made from a linear least-squares fit between these two quantities derived 
from those for galaxies for which both measurements were available. The conversions between the fully corrected line 
widths are as follows: 

For the Bottinelli et al. (1990) sample: W50 = 0.98 * W20 - 19.21 [N = 3325] 

For the Theureau et al. (1998) sample: W50 = 1.00 * W20 - 17.01 [N = 2055] 

Haynes et al. (1999) give a different measure of the line width denoted W21 and defined as the full width between 
the velocity channels at the 50% level of each horn. We can convert the Haynes et al. (1999)'s W21 to Bottinelli et al. 
(1990)'s W50 as: Wb0(Bot) = 1.02 * W21{Haynes) + 1.39 [N = 275], 

and Courteau (1997)'s Vmax and Botinelli et al. (1990)'s W50 are matched with: 

W50(Bot) = 0.99 * y max - 9.31 [N = 122] 

For 4 of the LSBs in the BdJOO sample we obtained V max values from de Blok, McGaugh, & Rubin (2001) and 
de Blok & Bosma (2002). A comparison between the Courteau (1997) V max values and the Bottinelli et al. (1990) 
IP50s for galaxies in common reveals little difference, and we assume here that Vmax values for the LSBs also map 
directly into W50. 

Multiple measurements of a galaxy from differnet samples vary typically by less than 10%, more than accurate 
enough for our purposes. In order to ensure accurate V ro t measurements, we restricted ourselves to galaxies with 
inclinations greater than 35°. Our final sample has 98 galaxies with reliable rotational velocity measurements. 
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Fig. Al. — Correlation between galaxy parameters: rotational velocity, V ro t (kpc), central surface brightness ((qj jj, total magnitude, 
AfH m ,Ki an< i disk scale length hu.K- The dotted line is at logio(V / r ot = 120 kms -1 ). Points with dotted error bars represent the galaxies 
from our sample. Points with solid error bars are the BdJOO sample. 



Fig. A2. — Correlation between galaxy parameters and morphological type. The dotted line is at logio(V ro t = 120 kms 1 ). 



Fig. A3. — Comparison of the stellar population model tracks for different initial mass functions; the Salpeter (1955) IMF (eq. 12*1, and 
the Chabrier (2003) IMF (eq. 0) for the GALAXEV models (left panel) and of the Salpeter versus the Kroupa (2001) (eq. |3j IMFs for 
the PEGASE models (right panel). The GALAXEV models are also plotted on the right panel for comparison between two different SSP 
models. 



Fig. A4. — Time evolution for the exponential (eq. ^J, upper panel, and Sandage (eq. 0, lower panel, star formation historie s (so lid 
curves). The dotted curve is a Sandage-style burst of star formation in which 10% of the total mass of stars are formed. See Figure IA7l for 
the effect of such a burst on the population model grids. 
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Fig. A5. — Comparison of the stellar population model tracks for exponential (eq. [BJ and Sandagc (cq.[5J star formation histories. 
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ta^iEjBjft^stT6etg^i,Stftg-iMi,gt^Iti6st#jit«see ilofecpietjj^mijpfnn^; ^m§Rft((fiigh&rfl^)-H m is the "modified" iT-band magnitude for 
the Courteau et al. sample, converted to A'-band with 2MASS H — K colors (see text for details) for direct comparison with the X-band 
data of the BdJOO sample. 



Fig. A7. — Comparison of the stellar population model tracks for a si mple exponential SFH (eq. |5J and one with a 4 Gyr old Sandage 
(eq. |SJ burst contributing 10% of the total mass (left plot, and see Figure |A4*1 and 50% of the mass (right plot). 

Fig. A9. — Near-IR-optical color— color plots for the Courteau et al. sample separated by morphological type. The galaxy centers are 
indicated by solid symbols: circles for Type-I galaxies, triangles for Type-II, and asterisks for Transition galaxies. Open circles denote the 
half-light radius of the disk. The line types for the galaxy color profiles represent bar strength: solid for bar-less (A), dashed for mild 
bars (AB), and dot-dashed for strong bars (B). Average errors due to the global calibration, as well as sky uncertainties for the innermost 
and outermost points are shown as crosses. Over-plotted are the Bruzual & Chario t (20 03; GALAXEV) stellar population models for an 
exponential SFH at different mctallicities. Lines of constant average age (A) (see eq llOl are the dashed, roughly vertical lines, while lines 
of constant metallicity are solid and nearly horizontal. In the upper left corner of each panel, triplex and foreground screen dust models 
are plotted. For the triplex model, the Milky Way dust extinction curve and albedo from Gordon, Calzetti, & Witt (1997) was adopted. 
The galaxy center and disk half-light radius are denoted by solid and open circles respectively. 
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Fig. A10. — Near-IR-optica l col or-color plots for the Courteau et al. sample separated by rotational velocity, V TO t (kms ). Symbols 
and line- types are as in Figure lAHl 



Fig. All. — Near-IR-optical color-color plots separated by rotational velocity, V IQ t (kms 1 ), for the BdJOO sample. 
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Fig. A12. — Near-IR- optic al color-color plots for the Courteau et al. sample separated by H-band magnitude, Mjj. Symbols and 
line- types are as in Figure IA9I 



Fig. A13. — Near-IR-optical color— c olor plots for the Courteau et al. sample separated by H-band central surface brightness, fiQ jj. 
Symbols and line- types are as in Figure lA~9l 



Fig. A14. — Average age, (A) (left), and metallicity, \og\o(Z / Zq) (right) as a function of local H m and i-sT-band surface brightness, where 
Hm is the "modified" H-band SB for the Courteau et al. sample, converted to i^-band with 2MASS H — K colors (see text for details) for 
direct comparison with the i-f-band data of the BdJOO sample. The H m data is distinguished with dotted error bar lines. The different 
symbols represent total Hm, i^-band magnitude bins for each galaxy. The dotted lines in the metallicity plot (right) denote the limits of 
the population model grids. Points at or near these limits for metallicity, and those at 12.9 Gyr for the average age should be interpreted 
with caution. 
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Fig. A15. — "Effective" average age (average of 0.5— l.5h bin), (A) e re as a function of central surface brightness (left) and total H m , it-band 
galaxy magnitude (right) for all 158 galaxies. See caption of Figure IA14I for the definition of H m . 

Fig. A16. — "Effective" average age as a function of scale length (left) for all 158 galaxies and rotational velocity (right), V TO t for all 98 
galaxies for which we have reliable V TO t measurements. 
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Fig. A17. — "Effective" metallicity logio(^ e ff /Zq) as a func tion of central surface brightness (left) and total if m , if-band galaxy 
magnitude (right) for all 158 galaxies. See caption of Figure IA14I for the definition of H m . 



Fig. A18. — "Effective" metallicity, logio(Z c ff /Zq), as a function of scale length (left) for all 158 fitted galaxies, and rotational velocity 
(right), V TO t for all 98 galaxies for which we have reliable V TO t measurements. 
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Fig. A19. — Age as a function of radius for all galaxies. The left plot gives the radius in terms of the measure scale length hu jf while 
the right plot gives the radius in physical units (kpc). Dotted lines represent our sample whereas solid lines are the BdJOO sample. 



Fig. A20. — Same as Figure lA"l9l but for mctallicity. 
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Fig. A21. — Age gradient as a function of central surface brightness, Mo,ir m ,if f° r al l galaxies. The gradients are fit out to a different 
number of radial bins: 0-1.5 hu,K (top left), 0-2.5 hu x (top right), 0-3.5 hfi,K (bottom left), and 0-4.5 hjjj^ (bottom right). 
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Fig. A22. — Age gradient as a function of Hm,K-ba,nd central surface brightness, ^to,ff m K (left) and total magnitude, Mjj m k (right) for 
all galaxies. The fit ranges represent "inner" gradients (left panels) fit out to 1.5 H ,K-band disk scale lengths, hjj,K^ and "outer" gradients 
(right panels) fit from 1.5 to 3.5 hu,K- The upper and lower plots show gradients as a function of scale length and kpc respectively. 



Fig. A23. — Same as Figure lA22l except for disk scale lengths, log(/i£f k) (left) , and rotational velocity, V ro t (right). 
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Fig. A24. — Same as Figure lA22l except for metallicity gradients. 



Fig. A26. — Metallicity gradient as a function of age gradient for all 158 galaxies. 



Fig. A27. — Effective age (left) and metallicity (right) as a function of morphological type for all 158 galaxies. 
Fig. A25. — Same as Figure lA23l except for metallicity gradients. 
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Fig. A28. — Age gradient (left) and metallicity gradient (right) as a function of morphological type for all 158 galaxies. The upper panels 
plot radius in terms of the measured scale length hn,K while the bottom panels plot radius in physical units (kpc). 



Fig. A29. — Age gradient (left) and metallicity gradient (right) histograms as a function of barredness for all 158 galaxies. The upper 
panels plot the radius in terms of the measured scale length hn,K while the bottom panels plot the radius in physical units (kpc). The 
fit ranges represent "inner" gradients (left panels) fit out to 1.5 hif,K, an d "outer" gradients (right panels) fit from 1.5 to 3.5 hn,K- The 
vertical dotted lines are located at zero gradient for reference. 



